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CLUMPY STREAMS FROM CLUMPY HALOS: 
DETECTING MISSING SATELLITES WITH COLD STELLAR STRUCTURES 
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ABSTRACT 

Dynamically cold stellar streams are ideal probes of the gravitational field of the Milky Way. This 
paper re-examines the question of how such streams might be used to test for the presence of "miss- 
ing satellites" — the many thousands of dark-matter subhalos with masses 10 5 — 10 7 M Q which are 
seen to orbit within Galactic-scale dark-matter halos in simulations of structure formation in ACDM 
cosmologies. Analytical estimates of the frequency and energy scales of stream encounters indicate 
that these missing satellites should have a negligible effect on hot debris structures, such as the tails 
from the Sagittarius dwarf galaxy. However, long cold streams, such as the structure known as GD-1 
or those from the globular cluster Palomar 5 (Pal 5) are expected to suffer many tens of direct im- 
pacts from missing satellites during their lifetimes. Numerical experiments confirm that these impacts 
create gaps in the debris' orbital energy distribution, which will evolve into degree- and sub-degree- 
scale fluctuations in surface density over the age of the debris. Maps of Pal 5's own stream contain 
surface density fluctuations on these scales. The presence and frequency of these inhomogeneities sug- 
gests the existence of a population of missing satellites in numbers predicted in the standard ACDM 
cosmologies. 

Subject headings: cosmology: theory - dark-matter - Galaxy: halo - Galaxy: kinematics and dynamics 
- Galaxy: structure 



1. INTRODUCTION 

Standard ACDM models of the Universe allow us to ex- 
plain structure formation on large scales. However, they 
predict an order of magnitude more dark-matter subha- 
los within the halos of typical galaxies than the num- 
ber of known satelli te galaxies orbiting the Milky Way 
(IKlypin et all 1199ft IMoore et al.l 1199ft IDiemand et~aLl 
120071: ISnringel et all 12008^ 7" Recent, large-area stellar 
surveys have discovered dozens of new satellite galaxies, 
mos t notably using the S l oan Digital Sky Survey (SPSS , 
e.g. iWillman et all [20051 IBelokurov et all l2006al 12007 ' 



Zucker et al. 20061: llrwin et al.l l2007t IKoposov et^il 



2007 1 IWalsh et all [2007D but "the number discrepancy 



between simulated dark-matter subhalos and observed 
satellite populations is still significant. This discrep- 
ancy can partially be explained by accounting for the 
the incomplete sky-coverage of SDSS and the distance- 
dependent limit on t his survey's sensitivity to low-surface 
brigh tness objects (K oposov et al.l 120081 : iTollerud et al.1 
2008). Indeed, models which take this into account and 
consider diffuse, (i.e. undetectable) satellit e galaxies can 
recon cile the number counts for subhalos (Bullock et al. 
2010). However, when they impose the suppression of 
stellar populations in low mass subhalos (which have 
masses below 5 x 10 8 M Q ) the number of undetected 
galaxies significantly declines and the prediction of nu- 
merous purely dark-matter subhalos less massive than 
5 x 10 8 M Q remains. 

There could be a genuine absence of "missing satel- 
lites" in the inner halo due to destruction by disk 



shocks , as illustrated in the calculations D'Onghi a et al.1 
(2010). However, note that these analytic descrip- 
tions of disk shocking based on the energy criterion 
are known t o overestimate disr uption rates of subhalos 
significantlv (|Goerdt et al.l[2007l ). Once these destructive 
effects are accurately accounted for, proof of the exis- 
tence (or lack) of these "missing satellites" could provide 
an important constraint on the nature of dark matter, 
which sets the minim um scale for the for mation of dark- 
matter subhalos (e.g. iHooper et a l. 2007). 

Along with the discovery of new satellite galaxies, 
SDSS has also uncovered a multitude of stellar struc- 
tures in the Milky Way halo from disrupting glob- 
ular clusters or satellite galaxies. In many cases, 
the debris is dynamically cold and distributed nar- 
rowly i n space (lOdenkirchen et al.ll200ll IBelokurov et all 
2006 k lLauchner et all 120061 : iGrillmair fc Johnson! 120061 : 
Grillmair fc Dionatosl2006btlGrillmairl2006Ll2009D . Such 
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cold stellar streams should be sensitive probes of the 
gravitational potential. On global scales, they can be 
used to constrain the radial profile, shape and orienta- 
tion of the Milky Way's triaxial dark-matter halo (e.g 
Johns ton et al.l 1199ft llbata et all 120021: Uohnston et al 
20051: iBinnevI l2008t lEvrd 120101 : fKoposov et all 1201 
Law fc Maiewskil I2010D . The presence of dark-matter 
subhalos would add asymmetries to the global potential 
over a range of smaller scales which will perturb these 
cold streams or even destroy them. Hence, if the missing 
satellites do exist they will add random uncertainties to 
any stellar-dynamical assessment of the global potential. 

Cravitational lensing has been suggested to be 
a u seful tool t o probe the presence of subha- 
los dChibal [20021: IMetcalf fc Zhaol 120021 iChen et all 



20031: iMoustakas fc MetcalJ 120031: IMetcalf et all 120041: 



Keeton fc Moustakas 2009; Rich m et al.ll200ft Ku et all 
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2009). These investigations conclude that flux ratio 
anomalies in lensed images or lensing time delays could 
be caused by dark-matter subhalos, though the con- 
straints are limited by our knowledge of the spatial dis- 
tribution of subhalos. However, this method is only ap- 
plicable to the most massive and most centrally concen- 
trated dark-matter halos, and not to galaxies like the 
Milky Way more generally. 

The effect of dark-matter subhalos on ste llar streams 
has be en explored in several previous studies. Ilbata et al.l 
(2002) showed that debris from the destruction of a 
10 6 M© globular cluster should be affected by heating 
due to repeated close encounters of subhalos and con- 
cluded that this effect could be detectable w i th fu - 
ture astrometric surveys. Moreover, iQuinn et~aT1 (2008) 
found that the inhomogeneities seen in Pal 5's tidal tails 
could not be accounted for in simulations evolved in 
a smooth potential. For the streams of larger satel- 
lites like the Sagittar ius dwarf galaxy (hereafter Sgr), 
Uohnston et al.l (|2002[ ) found that although stars in the 
debris are scattered by encounters with dark matter sub- 
halos, the thickness of the current Sgr stream could be 
explai ned as being due to the Large Magellanic Cloud 
alone. iSieg al-Gas kins fc Valluril (|2008t ) tested the addi- 
tional influence of different host potentials on debris from 
satellite galaxies and pointed out that while subhalos can 
shift the positions of streams and cause clumpy struc- 
tures, the shape of the halo potential and orbital path 
can have an ov erall comparable effect. Most recently, 
ICarlberel (|2009f) modeled a simplified stream on a circu- 
lar orbit and concluded that dynamically old (> 3Gyr) 
stellar streams cannot survive in the presence of subhalos 
with the masses and numbers predicted by ACDM. 

These previous works point to stellar streams as per- 
haps the most powerful way to find the missing satellites. 
However, none of these investigations separated the effect 
on streams of the known (and therefore uninteresting!) 
satellites with masses > 10 8 M© (jBovill fe Ricottil 120091: 
lBullockll2010L and references therein) from those that are 
"missing" (the pure dark- matter subhalos). In this study 
we construct a framework for understanding stream in- 
homogeneities by first isolating and dissecting the char- 
acteristics of disturbances caused by dark-matter subha- 
los alone. In contrast to previous work, which looked 
at the overall response of cold streams to the complete 
ACDM subhalo mass spectrum, we look at the expected 
frequency, influence and characteristic observable signa- 
tures of subhalos in each mass decade separately. We also 
contrast the response of different streams to the same 
masses, from ribbons such as Pal 5 to the giant stream 
from Sgr. Our twin aims are: (i) to understand with 
which streams we are most likely to be able to conclu- 
sively prove the existence or absence of missing satellites: 
and (ii) to learn how signatures of missing satellites that 
are apparent in streams might be interpreted. 

It should be noted that the discovery of very cold 
streams from globular clusters has inspired discussions of 
how the intrinsic properties of stellar streams themselves 
could cause inhomo geneit ies in their density distribu- 
tions (jKiipper et al.1 120081 l2010t iQuillen fc Comparettal 
120101) and these self-induced fluctuations could confuse 
the conclusive association of observed disturbances with 
dark-matter subhalo interactions. Our own work is 



also motivated by these current observations which con- 
tain tantalizing suggestions of non-uniformity in some 
cold stellar streams (e.g. the structure known as 
GD-1 and those from the globular cluster Pal 5, s ee 
l Odenkirchen et al.l 120031 : iGrillmair fc Dionatosl l2006al lbl: 
iKoposov et al.ll2010( ) , as well as the prospect of the den- 
sity of these streams being mapped more extensively (and 
accurately) in space and velocity with observations in the 
near future. Such observations could potentially distin- 
guish between the effect of subhalo encounters on differ- 
ent mass scales as well differentiate these signatures from 
non- uniformities due to intrinsic stream dynamics. 

We first review our understanding of the properties of 
dark-matter structures and stellar stream evolution in 
smooth potentials in <j2j We use this understanding to 
make order-of-magnitude estimates for the frequency and 
effect of encounters of stellar streams with structures of 
different masses in Sj3l We then go on to illustrate these 
expectations with numerical experiments in Sf5] as well as 
discuss the observational signatures of these encounters 
in <|5l We summarize our conclusions in §6. 

2. BACKGROUND AND METHODS 

The aim of this paper is to characterize the effects that 
dark matter subhalos orbiting around the Milky Way can 
have on debris from satellite disruption. In order to iso- 
late the influence of subhalos from other factors , neither 
the debris distribution nor the Milky Way is modelled 
self-consistently, and a simplified form for the mass dis- 
tribution is assumed in both cases. In particular, in our 
numerical experiments, the debris is represented by 3,000 
test particles which orbit, along with a varying number 
of subhalos, in a smooth and spherical Milky Way halo. 
The orbits of the test particles and subhalos are inte- 
grated using the leap-frog method with a time-step of 
0.5Myr. The test particles respond to the gravity of both 
the Milky Way halo and subhalos, but the subhalos do 
not interact with each other. 

Note that our numerical approach misses several ef- 
fects which could themselves contribute to non-uniform 
appearance of streams. These include multiple or contin- 
uous mass-loss, other asymmetries in the Galactic poten- 
tial and self-gravity of th e sta r streams. We will discuss 
each of these further in § 15.31 

The dark-matter distributions (both parent and sub- 
halos) are chos en to mimic the end point of the Via 
Lactea II (VLII. lDiemand et al. 2008) N-body simulation 
of structure formation on Galactic scales (see $2.1]\ , while 
the distribution of tidal debris (and test particle orbits in 
our numerical experiments) are chosen to match expec- 
tations from N-body simulations of satellite destruction 
(see $L2]). 

2.1. Dark Matter Halo Properties 

The parent dark-matter distribution is a Milky- Way- 
like halo r epresented by a Nav arro- Frenk- White (NFW) 
potential ([Navarro et al.l[l996h : 

p{r) = S c 

Pent (r/r s )(l + r/r s ) 2 

with parameters chosen to match the VLII N-body sim- 
ulation (M MW = 1-77 x 10 1 2 M w ,fi s = 24.6kpc, R V1I = 
389kpc, iDiemand et"aI1l2008D . 
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Figure 1. Number, scale radii (r s ) and tidal radii (r t id a i) of sub- 
halos in VLII as a function of subhalo mass (M su b). Gray points 
are for all subhalos and black are for those within 25 kpc at the 
present day. 

We also assume the NFW form for the subhalos and 
directly take all physical properties (i.e. the masses, 
tidal and scale radii, see Figure [J) and orbits of subhalos 
from t he publicly available analysis of the VLII simu- 
lation (|Diemand et al.lfeOOSf ). In practice, cosmological 
simulations like VLII have finite resolution, so the mass 
function is not complete below 4 x 10 6 M Q . We com- 
pensate for this resolution limit by duplicating subhalos 
within the mass range of 1O 5 M0 - 10 7 M Q to maintain 
the power-law mass spectrum with a power law index 
-1 set by larger masses (see upper panel of Figure [T]). 
As a result, each mass decade has 10 times more sub- 
halos than the mass decade above. The orbits of the 
duplicated subhalos are set by rotating the position and 
velocity vectors of the original subhalos by random an- 
gles so as to preserve spherical symmetry and velocity 
isotropy of the subhalo population. This method retains 
the potential and kinetic energy of subhalos and their 
radial distribution. 

2.2. Tidal Stream Properties 

The general principle of debris evolution can be sim- 
ply described: stars being torn from a satellite initially 
share a single orbital phase (or angle), but have a small 
range in orbital properties (or actions); this range in or- 
bital properties corresponds to a range in orbital peri- 
ods (or frequencies); and these period differences lead to 
gradual sp reading of the d e bris along the satellite 's orig- 
inal orbit () Johnston! fl998t iHelmi fe White! 1 19991 ) . Note 
that this description of evolution as simply phase-mixing 
can break down for larger mass satellites (more than 
0.01% of the mass of the parent) whose own gravity has 
been shown to in fluence the morphology of the debris 
(jChoi et all 12001 . This effect is ignored in our work 
since we are mostly interested in lower mass satellites 
and in differences in debris distribution at smaller scales. 

2.2.1. Scales in Debris 

Uohnstonl (|1998l ) and lJohnston et all (|2001ft found that 
the distribution of debris particles observed in N-body 



simulations of satellite disruption could be described in 
terms of the dimensionless tidal scale 



1/3 



(2) 



where m is the mass of a satellite and M p is the mass of 
a host halo enclosed by the pericenter, R p , of the satel- 
lite orbit. The tidal radius where the internal and tidal 
forces are in equilibrium and mass loss occurs scales as 



ftide 



sR p = 



1/3 



(3) 



The orbits of debris are offset in energy, E, and angular 
momentum J from the satellite's own orbital properties. 
The tidal scales can be used to derive order-of-magnitude 
estimates for the characteristics sizes and ranges of these 
offsets, AE ~ e and A J ~ j, where 



rtide 



and 



dR 



= sJ. 



GMp 
R„ 



(4) 



(5) 



Given these characteristic ranges, the angular length 
of streams of debris as viewed from the Galactic center 
is expected to grow by of order 



A* = e 



2ir dTy 



Tq, dE 



(6) 



each orbit, where J c i rc is the angular momentum of a cir- 
cular orbit of energy E and azimuthal time-period Tq,(E) 
at radius R c i rc (E). The secondary dependance of orbital 
time periods on angular momentum has been ignored. 
The angular width is initially of order 

w = s. (7) 

and is expected to gro w with time at a rate dependent 
on the parent potential (IHelmi fe Whitd[l999l ). For near- 
spherical potentials the rate is sufficiently small that the 
approximation w ~ s remains useful for many orbits 
(|Johnston et al.H200lD . 

2.3. The case of Palomar 5 

Tidal streams associated with the globular clus- 
ter Palomar 5 we r e dis c overed using SDSS data by 
lOdenkirchen et all (|200lL I2003D stretching several de- 
grees away either side of the center of the cluster, and 
have now been mapped to a total extent of 22 degrees 
(jGrillmair fc Dionatosl[2006a|) . They were the first de- 
bris to be mapped from such a small object and remain 
a primary example of a thin, cold stream. 

The first set of rows of Table Q] list the observed char- 
acteristics of Pal 5's stream. Since Pal 5's orbit is not 
much further from the Galactic center than the Sun, 
these are roughly translated into a Galactocentric view 
in the second set of rows, assuming the Sun is at 8 
kpc from the Galactic center — for example, a lower 
limit on the Galactocentric angular length is taken to 

be I = (-Ddebris/-RdebrisKobs where Z?dobris and -Rdcbris 
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Object 


Pal 5 


Sgr 


Heliocentric view 


£>dcbris (kpc) 


~ 23.5 


8-80 


i"obs (degrees) 


0.5 


10 


'obs (degrees) 


22 


>360 


Galactocentric view 


-Rdebris (kpc) 


18 


8-80 


/ (degrees) 


> 29 


> 360 


assumed properties 


R p (kpc) 


7.5 


15.0 


i?a (kpc) 


19.2 


60.0 


Tr (Gyrs) 


0.33 


0.89 


T* (Gyrs) 


0.55 


1.35 


Rene (kpc) 


13.7 


40.1 


m S at(M ) 


10 4 


5 x 10 8 


age (Gyrs) 


8.44 


1.9 


derived scales 


s 


0.007 


0.17 


e (km/s) 2 


123 


4572 


A>1/ (degrees) 


1 


42 


w (degrees) 


0.4 


10 



Table 1 

Properties of the Pal 5 and Sgr streams. 



are the heliocentric and Galactocentric distances respec- 
tively. The observed properties can be broadly repro- 
duced with an N-body simulation of a hot, spherically 
symmetric system disrupting along an orbit with charac- 
teristics listed in the third set of rows ("assumed prop- 
erties") of Table [TJ The simulation adopted a 10,000- 
particle Plummer model with mass 10 4 M Q and scale 
lengt h 7.5 pc for Pal 5, and calc ulated its self-gravity 
using Hcrnauis t fc Ostr ikcr (1991) "self-consistent-ficld" 
code. The model was allowed to evolve in the NFW 
parent potential described in Section 12.11 for 8.44 Gyrs. 
Very similar masses an d orbits were found w ith the more 
extensive modeling by D ehnen et al.l (|2004f) . 

The left-hand panels of Figure [2] summarize the results 
of our simulation . The top left-hand panel shows the 
energy- angular momentum distribution of debris parti- 
cle, with axes scaled by the factors given by equations 
id} and iJSJ), and listed in the bottom set of rows in Ta- 
ble [TJ Note that debris particles are systematically offset 
from those still-bound to the satellite, and distributed 
in the range ±3e and ±3j around the satellite's own or- 
bital energy /angular- momentum — a result that is found 
to be largely inde pendent of sate llite mass, scale, pro- 
file and orbit (e.g. lJohnstodlT998l ). Note that the range 
in scaled-angular-momenta does depend (mildly) on the 
eccentricity of the orbit, with debris on more eccentric 
orbits exploring larger ranges. 

The bottom left-hand panel of Figure [2] shows the final 
position of the particles projected onto the plane of the 
sky, as viewed from a Solar position in our model. From 
our assumed mass and orbit, our analytic estimate (i.e. 
equation [7]) suggests an angular width w ~ 0.4 degrees 
for the debris. While this is of the same order as both 
the observed and simulated width, this prediction is rel- 
atively insensitive to the assumed mass for the satellite, 
so the agreement should only be seen as a weak confir- 
mation that the mass of the satellite is the correct order 
of magnitude. 

The estimated angular length of a stream of age t 
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Figure 2. Different representations of the globular cluster, Pal 5. 
N-body (left hand panels) and test particle streams (right hand 
panels) shown in energy /angular-momentum space (top panels) 
and in the sky (bottom panels). The dashed lines in the top pan- 
els show the mean orbital properties adopted for the test-particle 
streams. In the right-hand panels, the gray dots repeat the N-body 
results from the left and the black dots are the test particles. The 
dashed lines in the bottom panels represent the orbital path and 
the arrow shows the direction of the moving stream. 



viewed from the Galactic center is 

Rcivi 



I 



t 



Rd 



cbris 



A*, 



(8) 



where A^f is the expected angular growth per orbit cal- 
culated from equation (JSJ) for the given mass and orbit. 
The additional scaling, dependent on the radial phase 
at which the debris is observed (i.e. distance of debris 
from the Galactic center i?debris), is included to account 
for the effect of debris density increasing/decreasing as 

the angular speeds decrease/increase (oc (l/i?dcbris) 2 in 
a spherical potential from conservation of angular mo- 
mentum. The extra factor of 4 reflects the range that 
characterizes the width of the energy distribution seen 
in the top left-hand panel of Figure [2] (i.e. ±2e). The 
prediction for the angular length of the stream, as ob- 
served from the Solar position, is given by 



^obs — 



cbris 



D, 



debr 



(9) 



For the parameters adopted in our simulation of Pal 5 
we find lobs ~ 30 degrees, which agrees with the length 
seen in the lower- left panel of Figure [2] 

2.3.1. Initial Conditions for Numerical Experiments 

This paper employs a simplified representation of N- 
body results that characterizes the evolution of tidal de- 
bris described above using test particles integrated in 
the combined potential of the parent and orbiting lumps, 
but ignoring the self-gravity of the satellite. The parti- 
cles are given initial positions and velocities slightly offset 
from the satellite's own, that are chosen to reproduce the 
scaled energy and angular momentum distributions seen 
in the top left-hand panel of Figure [2j Specifically, the 
particles are initially positioned uniformly distributed 
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along two lines pointing in the direction of the Galac- 
tic center and centered on the satellite at the apocenter 
of its orbit at Galactocentric radius R = R a and tangen- 
tial velocity Utan- Particles at position AR along the bars 
relative to the satellite are assigned a tangential velocity 
Utan + Av, where these offsets in position and velocity 
satisfy the equations: 



/ d$> \ 

AE = ae = v tan Av + AR — 

\dRJ R=R a 

AJ = /3 j = VtanAi? + RAv. 



(10) 



Here a is restricted to the range ±3 and P = 2.5a/3 in 
order to mimic range in the orbital properties seen in the 
simulations In addition, the particles are assigned small 
random velocities in the x, y and z directions around 
the mean u ta n + Av chosen from a gaussian with width 
a = 0.6swtan,apo to reproduce the spread in angular mo- 
mentum at a given energy. 

The right-hand panels of Figure [5] illustrate the success 
of this method by comparing energy/angular- momentum 
distributions (top panels) and final positions (bottom 
panels) for our N-body (gray particles) and test parti- 
cle (black points) model of debris evolution for Pal 5 
(assumed properties in third set of rows of Table [T]). 

3. RESULT I: ANALYTIC ESTIMATES 

In this section we use simple analytic arguments to de- 
rive order-of-magnitude estimates for the frequency of en- 
counters of streams with varying mass dark-matter sub- 
structures (333J and the characteristic scales of these 
subhalos' influence on the debris f £|3.2[) . We first illus- 
trate the implications of these estimates for the debris 
from Pal 5 and then go on to discuss the application to 
other streams (JO]). 

3.1. Frequency of Encounters 

Suppose a (perfectly tubular!) stellar stream of angu- 
lar length I and angular width w is orbiting at an average 
distance R clrc from the Galactic center through a sea of 
dark-matter subhalos of size r t idai and number density 
^subhaio, The rate of encounters that directly impact the 
stream with relative speeds in the range (v,v + dv) can 
be estimated as: 



dN r 



dt 



volume encounter rate x number density x P re \(v)dv 

■ (stream surface area) x v x n su bhaio x Prc\{v)dv 

■ (i? circ I x 27r& max ) x n sub haio x v P rc i{v)dv. (11) 
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where P le i(v) is the probability distribution of relative 
speeds. The maximum impact parameter for direct en- 
counters is taken from the condition that the stream and 
subhalo overlap physically: 



Figure 3. Encounter frequencies. The upper panel plots the es- 
timated number of direct encounters for Pal 5 (solid lines) and Sgr 
(gray lines) debris during the debris lifetime and for subhalos of 
different masses. Each line presents the case for different subhalo 
mass decades as noted on the lines. The middle panel shows the 
distribution of encounter velocities. The bottom panel shows an 
estimate for the fraction of encounters of fixed energy that are due 
to slow encounters with small masses (see text). 



MAX KR circ /2, nidai]. 



(12) 



Assuming that the local velocity distribution of the sub- 
halos is an isotropic Maxwellian with dispersion cr, we 
can approximately represent the relative distribution of 
encounter speeds along a typical orbit by the relative 
speed distribution between the subhalos 



P Te i(v)dv 



1 



(V2^r) 3 



v 2 exp 



(see lBinnev fc Tremainj|2008l equation (8.46)). This im- 
plies a typical encounter speed of 



„typ — 



/ °°^ 2 P rc i(^)^ 
f£° vP ie \(v)dv 



= ZsfKo/2. 



(14) 



Replacing I in equation (TTTj) with equation ([5]) (evalu- 
ated for -Rdebris = -Scire) and integrating over speeds and 
time, we find the number of encounters with encounter 
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speeds less than w cnc that this stream experiences over 
its lifetime t: 

N cnc (v < v cnc )=8^/TrR cilc b metK at n sub haio X A* (-^- 



l-fl + ^)ex P 



Aa 2 



.(15) 



The solid lines in the top panel of Figure [3J show 
N enc (y < v enc ) for our Pal-5-like stream (i.e. evaluat- 
ing equation [15] with parameters adopted from Table 
[IJ as a function u cnc for different decades in subhalo 
masses and assuming a characteristic a of 120 km/s. 
(Note that the circular velocity at radius 25 kpc in a 
spherical NFW halo with parameters chosen to match 
the parent halo in VLII is 178 km/s, which suggests 
a dispersion of 178/-\/2=126 km/s.) The subhalos in 
each group are assumed to have tidal radii similar to 
those seen in VL II for subhalos within 25kpc (i.e. 
?"tidai = 0.21, 0.58, 1.6, 4.5, 12. 7kpc for masses in ranges 
10 5 - 10 6 /10 6 - 10 7 /10 7 - 10 8 /10 8 - 10 9 /10 9 - 1O 1O M 
respectively — see Figure [TJ . For n S ubhaio we adopt the 
mean number density within 25 kpc given by assuming 
there to be 42 subhalos in the rang e 10 6 - 1O 7 M within 
this volume (as seen in VLII, see Figure [J), and a factor 
of 10 more/less in each subhalo mass decade below/above 
(to mimic the mass spectrum seen in VLII). (This is ex- 
pected to be an upper limit on the number of encounters 
since the adopted age produces a stream slightly longer 
than the observed stream, and the mass adopted for our 
Pal 5 model was fairly low.) Inserting these numbers for 
the 10 6 — 10 7 M Q mass range explicitly into equation [15] 
and letting v cnc — !• oo gives the total number of direct 
encounters: 



total 



(M sub = 10 6 - 1O 7 M ) = 20 ( -^-) 



( a 



Vl3.7kpcy 
t " 2 
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Overall, our calculations indicates that Pal 5's stream 
will have suffered hundreds of direct encounters with 
subhalos less massive than 10 5 M Q , ~ 70 with subhalos 
masses in the range 10 5 — 10 6 Mq, ~ 20 with subhalos 
of 10 6 - 10 7 M. o , ~ 5 with subhalos of 10 7 - 1O 8 M , 
and a few with subhalos of 10 s — 10 9 M©. Subhalos in 
the higher mass bins are very unlikely to directly hit 
the stream, although they will influence it through more 
distance encounters. In addition, note that while there 
are sufficient numbers of smaller subhalos to fully explore 
the relative velocity distribution, the few encounters with 
large subhalos are likely to take place close to the typical 
encounter speed around the peak of the velocity distri- 
bution shown in the second panel. 

3.2. Effect of Encounters 

In order to develop some understanding of the effect 
of subhalos on streams, consider the idealized case of a 
subhalo encountering a perfectly straight stellar stream 



aligned with the x-axis, with impact parameter b along 
the y-axis and relative velocity v cnc = (v x ,0,v z ). Over- 
all, the change in speed of the stars perpendicular to the 
stream is zero by symmetry. Assuming the change in the 
relative speed for stars in the stream is negligible during 
the encounter ( Av x << v x ) and simply integrating the 
equations of motion we find 



Av x = 



a x dt 

GM suh (v x t + x) 

-oo [(v x t + x y + b 2 +v 2 z t^/- 

2GM suh x v 2 

Vencb 2 V 2 (x/b) 2 + 

^cnc 



■dt 



(17) 



where x is the position of the star in the stream at the 
moment of closest approach of the subhalo, relative to 
the impact point. If the space motion of the stars in the 
streams in the Galactic rest- frame is ^stream, then the 
energy change is 

AE = v strcam -Av+iAu 2 = v stIcam Au x +0(Av 2 x ) (18) 

For direct encounters with subhalos of mass M su b : 

^stream 



AE(b = 0) = 2- 



(19) 



3.2.1. Lowest Mass Subhalos 



Figure [3J suggests Pal 5 experiences thousands of en- 
counters with subhalos less massive than 10 5 M Q during 
its lifetime. In this many, weak encounter regime the 
first term in equation (|18|) should average out to zero 
and leave the second term to heat the stellar stream. 
Despite their large number, since the heating term is 
second order, these encounters would have negligible ef- 
fect compared to the few encounters with larger subhalos. 
For example, while there would be a factor of ten more 
encounters with M su b ~ 10 4 — 10 5 M© than in the mass 
decade above, the net energy change due to each of those 
encounters would be a factor of 100 smaller, so the indi- 
vidual encounters with larger lumps would be generally 
be expected to have a more significant effect. 

Nevertheless, equation ([T9"l) indicates that the effect of 
a lower mass subhalo can be comparable to that of a 
subhalo in the mass decade above if the encounter with 
the lower mass is ~ 10 times slower. Indeed, the large 
number of encounters with low mass halos ensures that 
a small number will come from the low-end stream of 
the relative velocity distribution (see middle panel of 
Figure [3j and the effect of these few, slow encounters 
would not be expected to average out. The lower panel 
of Figure [3J compares the frequency of these slow en- 
counters with low mass subhalos relative to ones of com- 
parable influence in the mass decade above by plotting 
the ratio of the number of encounters of low mass sub- 
halos to higher mass subhalos of comparable influence: 
[10iV enc (i; < v cnc /10)]/N enc (v < w cnc ). Since this ratio is 
never much greater than 1% we conclude that while rare, 
influential low-speed encounters with low mass subhalos 
can occur, they will be much less frequent than encoun- 
ters of comparable influence with higher-mass subhalos. 
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Figure 4. Estimate for the scaled energy imparted by a direct 
encounter at x = at different points along a stream. The left 
hand panel shows the influence of subhalos of mass 5 X 10 5 /5 X 
10 6 /5 X 10 7 /5 X 10 s /5 X 1O 9 M . The right hand panel plots the 
result for a subhalo of mass 5 X 10 7 Mq, traveling perpendicularly 
to the stream and with impact parameters b = 0/2/4/8/16r s . 



3.2.2. Intermediate Mass Subhalos 

Pal 5 should have experienced ~20 direct encounters 
with subhalos in the mass range 1O 6 -1O 7 M , whose in- 
fluence is unlikely to average to zero. In order to gain 
some intuition for their effect, the left hand panel of Fig- 
ure H] plots equation (|T9l) for values for the stream and 
encounter speeds averaged over their respective velocity 
distributions (vj[^ eain = 2-\/2/7rer and = 3y / 7TCT/2.) 
and M suh = 5 x 10 5 /5 x 10 6 /5 x 10 7 /5 x 10 8 /5 x 10 9 M Q 
(solid... dot-dot-dot-dash lines), as a function of posi- 
tion x along the stream. The extended nature of these 
subhalos was taken approximately into account by re- 
placing M su b in equation ([TT]) with the mass enclosed 
within radius r = x for an NFW model with scale r s 
(M su b(r)). (Note that the integral was solved analyti- 
cally and did not take this into account). The subha- 
los had r s chosen to match those within 25kpc in VLII 
(r s = 0.05, 0.12, 0.30, 0.76, 1.94kpc respectively). In or- 
der to show the behavior of AE for these different mass 
ranges in a single plot, the energy change was normal- 
ized by the characteristic energy for a direct encounter, 
calculated by substituting x — r s in equation (111 



AE, 



char 



GM sub ^ r p cam 8V2GM ! 



sub 



, typ 

u cnc 



3tt 



(20) 



Most noticeable in Figure [4] is the abrupt sign-change in 
AE at the impact point — stars upstream of the impact 
point are pushed to orbits of higher energy and stars 
downstream of the impact point are pushed in the op- 
posite direction to orbits of lower energy. This effect 
exaggerates the energy gradient already present along 
the stream (see £j2.2[) and creates a gap in energy. The 
fraction of the stream affected by the encounter depends 
on the size of the perturber — with the smallest subha- 
los influencing only a subset of stream particles directly 
around the impact point (solid lines) and the largest 
subhalos influencing the whole stream (triple-dot-dashed 
line) . 

How influential these encounters are on debris evolu- 
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Figure 5. Scales in an encounter relative to stream properties 
for subhalos of different masses. The top panel compares Ai? c i lar 
(symbols) with the intrinsic energy spread in Pal 5's (bold line in 
the middle) and Sgr's (faint line on the top) debris. The middle 
panel compares the length of Pal 5's stream with the tidal radius 
of the subhalos (Sgr's stream is far longer than the plot). The 
lower panel gives an upper limit on the angular scales of gaps due 
to encounters with Pal 5 (bold diamond symbols) and Sgr (faint 
triangle symbols) debris. 



tion also depends on stream characteristics. The top and 
middle panels of Figure [5] assess this by plotting ASchar 
for subhalos of different masses and the spatial scale over 
which this energy change is important (taken to be the 
tidal radius rtidai)- For comparison, the energy scale over 
which debris orbits are spread for Pal 5 (i.e. e, calculated 
from equation|4]) and the physical length of the stream it- 
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self are shown as bold lines. These plots suggest that each 
direct encounter of Pal 5's debris with subhalo masses in 
the range 10 6 — 10 7 M Q should create gaps in energy that 
are a significant fraction of the spread of energy already 
present in the debris (i.e. A£? c har ~ e) for of order 10% 
of stars in the stream (i.e. r t idai ~ Z/10). 

Over time these gaps in energy will cause the stream 
stars around the impact point to spread apart from each 
other faster than stars in other regions of the stream, 
leading to physical gaps in the stream density distribu- 
tion. We can estimate an upper limit on the angular 
scale of these physical gaps by replacing the stream en- 
ergy range 4e in equation © with Ai? c har and calculating 
Zobs for Pal 5's lifetime and orbit (assumed parameters 
in Table 1). The bold diamonds in the lower panel of 
Figure [5] show this estimate. 

Overall, we conclude that direct encounters of Pal 5's 
debris with subhalo masses in the range 10 6 — 10 7 M Q 
during its lifetime could lead to about 20 gaps in density 
of order several degrees across, while masses in the range 
10 5 — 10 6 M Q could leave many more on sub-degree scale. 
Few, if any, direct encounters are likely to occur for the 
Pal 5 stellar stream for masses greater than ~ 1O 7 M0 (as 
shown in Figure [3]) , but could potentially lead to larger 
gaps. 

3.2.3. Indirect Encounters with Largest Mass Subhalos 

The right hand panel of Figure 0] repeats the same 
plot as the left hand panel, but this time for a 10 7 M© 
subhalo traveling perpendicularly to the stream (i.e. 
v cn c — t>enc = v z) an< ^ passing it at impact parameters b 
of 0,2,4,8,16 r s . Note that for these more distant encoun- 
ters the impulsive approximation under which equation 
(fl7|) was derived breaks down and the estimate is likely 
to overestimate the size of a subhalos' influence. The re- 
sults are nevertheless included here as a broad guide to 
general behavior, but should be treated with some skep- 
ticism. The effect of more distant encounters can be as- 
sessed more realistically using numerical integration (see 

ED- 

As the impact parameter increases, both the size of 
the energy change decreases and the scale over which it 
is felt grows. The energy imparted to the stream changes 
more smoothly as a function of position along the stream 
than in the direct encounter case, with the most extreme 
changes apparent at the edges. This is in part because 
the impact parameter is now comparable to the length 
of the stream itself and it is only at the edges that the 
velocity change is aligned with the stream velocity to 
produce maximum energy change. The net result is that 
the stream is smoothly stretched in energy, rather than 
forming an abrupt gap. It is this type of interaction that 
dominates for the larger masses: Figure [3] tells us we 
should expect few such encounters over Pal 5's lifetime. 

Note that the same behavior can be seen for individual 
encounters for smaller masses, but in those cases there 
are sufficient numbers of weak encounters for the first- 
order terms in energy change to cancel out. So we con- 
clude that indirect encounters of smaller masses are not 
important for leaving observable signatures. 

3.3. Summary and Application to Other Streams 

In summary, our analytic estimates indicate that Pal 
5's stream should be sensitive to dark subhalos, with 



unique signatures from different subhalo masses. The 
dominant influence on streams like Pal 5's will be from 
scattering of stars by encounters with subhalos in the 
mass decades 10 5 — 10 7 M Q where many tens of encounters 
occur during the stream's lifetime. These events will lead 
to gaps in the stream with characteristics scales reflect- 
ing the mass of the perturbers. Low-speed encounters 
from smaller subhalos can create comparable gaps, but 
these events occur much less frequently and so should not 
confuse the signal. Heating by perturbers in the lower 
mass decades should also be negligible. Direct encoun- 
ters with larger masses are unlikely (though catastrophic 
if they occur), and indirect encounters will tend to stretch 
rather than chop up a stream. 

Our estimates also suggest that streams with Pal 5's 
scales are the most interesting. The thin lines and sym- 
bols in Figure|3]illustrate the results of repeating the cal- 
culations made above for Pal 5's stream for Sgr's debris. 
In comparison to Pal 5, these streams have a far larger 
cross section for encountering dark-matter subhalos, but 
also have a far larger spread in orbital properties in which 
to hide the signatures of encounters. This means that the 
influence of encounters with subhalos less massive than 
10 8 Mq are likely to be negligible as: (i) both the spatial 
and energy scales of the change are much smaller than 
the inherent stream scales (see first and second panels of 
figure[5]), and (ii) the encounters become frequent enough 
that first order changes may be expected to cancel out. 
Subhalos more massive than 1O 8 M0 could leave their im- 
print on Sgr debris. However, these dark-matter struc- 
tures are generally thought to contain stars and hence 
should be detectable by other means. 

We conclude that for the purposes of finding missing 
satellites, cold streams from low mass objects should be 
the most sensitive probes. 

4. RESULT II: NUMERICAL ILLUSTRATIONS 

This section presents the results of scattering test par- 
ticle streams by subhalos with numerical integration. 
We present our results in orbital (i.e. energy /angular- 
momentum) space, coordinate space, and observable 
space. The Sun is assumed to be at (X,Y,Z) = (-8kpc, 
0, 0). 

In the following subsections, we look first at individual 
close encounters ( §4. lj) . then multiple encounters for sub- 
halos in each mass decade f £14.21) . and finally examine the 
combined effect of all subhalos in the full ACDM mass 
spectrumC? 14.31) . 

4.1. Individual Encounters 

We first follow a single artificial encounter with M BU b = 
10 6 - 5 M Q , r s = 0.118kpc and r t idal = 0.590kpc passing 
directly through a Pal 5-like stream perpendicular to 
the orbital plane with a relative velocity (v x ,v y ,v z ) — 
(0, 0, — 200km/s) to the stream. Figure [5] illustrate the 
effect of this direct encounter on the stream (fiducial re- 
sults in each panel labelled with (f)). The orbital plane 
of the stream is in the X-Y plane. The left column shows 
the moment of closest approach of the subhalo encounter 
moving along the z-direction and passing perpendicularly 
through the stream. The middle and right columns show 
the particle positions in the X-Y and energy /angular- 
momentum planes 4.34 Gyr after the encounter. There is 
a clear energy gap around the impact point in the stream 
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Figure 6. The various effects of encounters on a Pal 5-like stream. Our fiducial encounter (particles labelled (f) in each panel) has 
M au i3 = 1O 6 ' 5 M0, v cnc = 200km/s, b = 0. The top/middle/bottom rows show the responses when the velocity /impact /mass are allowed 
to vary around this. The "+" in each panel represents the Pal 5 center, and the circles show the size (r^ai) of the encounters and where 
they pass through. Note that three encounters are all overlapped in the top-left panel. The streams in the first, second, and last columns 
are centered at the Pal 5 center, impact point, and mean energy/angular-momcntum respectively. Streams other than the fiducial model 
are shifted by -3 and -6 along y-axis from their original positions in the second and last columns. 



that forms immediately after closest approach: the par- 
ticles behind the impact points (with higher energy) are 
accelerated (gain energy) and the ones ahead the impact 
points (with lower energy) are decelerated (lose energy) 
as predicted in ^3.2.21 This energy gap results in a spa- 
tial gap along the stream which grows with time due to 
phase mixing as illustrated in the middle column. Hence, 
a stream influenced by many encounters should end up 
having a clumpy structure produced by numerous energy 
and spatial gaps. 

Having illustrated the nature of the gaps in energy, we 
now investigate how these response of the stream to var- 
ied encountering conditions. The first row of Figure |J5] 
illustrates the effect of encounters with different veloc- 
ities relative to the stream (v CDC = 100, 200, 400km/s), 
the second row shows the gaps resulting from encoun- 
ters passing with different impact parameters (6i mpact = 
0x,4x,8 x r s which correspond to 0, 0.470, 0.940kpc), 
and the last row shows various mass encounters (M su b = 
lO 5 ' 5 , 10 6 - 5 , 1O 7 5 M ). Overall, when the velocity of the 
encounter doubles the gap becomes half, while smaller 
impact parameters and larger masses result in larger en- 



ergy gaps. These scalings agree with the trends in the 
equations (|17l) . (TT51) . and (fTTj)) . For a stellar stream or- 
biting with many dark matter subhalos, irregular clumps 
are expected in the stream that correspond to various 
energy gaps created by different encounters at different 
times. 

Lastly, the bottom panels suggest that while a small 
subhalo has an effect on the local region of a stream 
generating a small energy gap, a close encounter with 
a subhalo larger than 10 7 ' 5 M Q can globally distort the 
entire energy-angular momentum scales, since the size 
and influence of the subhalos are larger than the entire 
stream. This will be discussed in more detail in § 14.21 

4.2. Multiple Encounters 

To examine the integrated influences of many sub- 
halos, we perform separate simulations with subhalos 
drawn from each decade of the ACDM mass spectrum. 
We first illustrate our results by contrasting simulations 
with the ~30,000 intermediate mass subhalos in the mass 
range 10 6 — 10 7 M© and ~30 large mass subhalos of 
10 9 - 10 10 M Q . Figure [7] shows the final spatial and 
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Figure 7. The final distribution of stream particles from simula- 
tions of encounters with all the subhalos (upper particle distribu- 
tions) and the 20% closest encounters (lower particle distributions) 
in the mass range 10 6 - 10 7 M Q (left panels) and 10 9 - 1O 1O M 
(right panels). The top row shows the stream in the sky and 
the bottom row shows the energy/angular-momentum distribution. 
The arrows illustrate the direction of the stream movement and the 
cross represents the center of host satellite Pal 5. The grey lines 
illustrate the orbits of satellites. The lower particles have been 
shifted for clarity by -10° and by -4 along the y-axis in the top and 
bottom panels respectively. 



energy/angular- momentum distribution of these cases. 
First consider the overall morphologies in Figure [7J Re- 
call that the calculation in < i3. II suggests there should be 
~ 10 direct encounters with subhalos in the lower mass 
range, which is in rough agreement with the visual im- 
pression of the energy/angular-momentum distribution 
in the lower left-hand panel. In contrast, the stream 
with the large subhalos (the right column of Figure [Jj 
does not have any gaps. Instead, the orbital phase of the 
stream with subhalos in the mass range 10 9 — 10 10 Mq 
is shifted and the entire energy scale is changed as sug- 
gested in § 14.11 All panels in Figure [7] are plotted relative 
to a central particle (marked with +). The scaled energy 
q and angular momentum (A J / sJ) values of the central 
particles changed by (-0.74, -0.53) and (9.40, 4.00) for 
the streams with 10 6 - 10 7 M Q and 10 9 - 10 10 M o sub- 
halos respectively. This phase and energy shift will not 
be observable since there is no way to probe the original 
phase and energy. 

In order to disentangle the effect of direct (close) en- 
counters from indirect (distant) ones, we select the 20% 
closest encounters based on the minimum distance be- 
tween the subhalos and the stream particles that occurs 
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Figure 8. Final distribution of stream particles from simulations 
with subhalo encounters of varying mass ranges. The streams with 
no subhalos, subhalos in the mass range 10 5 -10 6 M Q , 1O 6 -1O 7 M , 
10 7 - 1O 8 M , and 10 s - 10 9 M Q are shown in the sky (top), radial 
velocity (middle), and energy /angular- momentum space (bottom). 
The stream with no subhalos is at the original position and the 
others are shifted down by 10 and 5 for the top and bottom panel 
respectively. The grey crosses represent the center of the host satel- 
lite. 



throughout the simulations. We then re-run the same 
simulation with only these subhalos present. The results, 
presented in the lower particle distributions in each panel 
of Figure [7J are almost identical to the ones with all sub- 
halos in each mass range, which suggest that heating by 
the more distant encounters is negligible. 

Figure [8] repeats the plots in Figure [7j for simula- 
tions with subhalos in all mass decades. While the 
stream without any subhalos shows a smooth morphol- 
ogy in the sky and energy/angular-momentum space, the 
streams with subhalos in the mass range 10 5 — 10 6 M Q , 
10 6 - 10 7 M Q , and 10 7 - 10 8 M o show eye-catching sig- 
natures such as spatial, velocity, and energy gaps. As 
expected, the smaller subhalos cause many more smaller 
scale clumps in the stream than the large ones. In con- 
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trast, the stream with 10 8 — 10 9 M Q subhalos does not 
have clear signatures. These subhalos are too large to 
leave small scale variations but too small to shift the en- 
tire stream as in Figure[7l Rather, these subhalos distort 
the energy-angular momentum distribution and slightly 
change the overall shape of the stream. 

The gaps in Figure M due to 10 5 - 10 6 M Q subha- 
los are more commonly smaller than those due to the 
10 6 — 10 7 M Q objects. This indicates that each mass 
decade might be separately detected. While small en- 
counters with subhalos will generally result in smaller en- 
ergy gaps which correspond to smaller physical gap than 
encounters with large subhalos, note that in the middle 
of the leading part of the stream with 10 5 — 10 6 M© sub- 
halos, there is a large gap (where ARA ~ —5°) which is 
as big as the ones in the stream with 10 6 — 1O 7 M sub- 
halos. On further investigation, we found this large gap 
can be attributed to a very slow encounter (which caused 
a large gap in energy) that occurred relatively early in 
the simulations (which allowed the gap to grow). How- 
ever, we expect this to be a rare occurrence as described 
in ij3.ll To confirm this expectation, we ran four addi- 
tional simulations with the same conditions as the origi- 
nal 10 5 — 10 6 M Q run but different initial starting points. 
In none of these cases did such large gaps appear. 

We can also extrapolate from these numerical results 
to consider the effect of the smallest "missing satellites" 
with M su b < 10 5 Mq. The scales of the gaps due to 
individual encounters with these subhalos will be even 
smaller than those due the mass decade above (like the 
ones apparent in Figure [Sj> and would thus require a very 
high resolution map to detect. Making such a map is ob- 
servationally challenging since Pal 5's stream contains 
only a similar number of stars to the globular cluster 
itself (i.e. of order 1,000 stars), and at much lower sur- 
face density compared to the backg round field popula- 
tion. Moreover, as argued in § 13.2. H and also apparent in 
Figure [HI the averaged heating of many such encounters 
is negligible compared to responses to higher mass sub- 
halos. We conclude that it is unlikely that the smaller 
subhalos would leave an observable signature in Pal 5's 
stellar streams. 

4.3. Stellar Streams in KCDM 

To assess the integrated effect of subhalo encounters, 
we would ideally run simulations of streams evolving 
in the presence of dark-matter subhalos from the full 
ACDM spectrum. However, as outlined in ^4.2[ close en- 
counters play a decisive role in shaping streams. Hence, 
we limit our computational expense by selecting only the 
97,090 subhalos which enter within a Galactocentric ra- 
dius of 25kpc at least once during the 8.44Gyr fiducial 
simulation time from our full realization of over 300,000 
subhalos in the mass range 10 5 — 10 10 M Q . 

Figure UH summarizes the results of 12 such simulations 
with projections of the final particle positions, centered 
on the satellite body, onto the x-y plane. All the simu- 
lations were of the same stream, but started one radial 
orbital period apart in order to explore a variety of en- 
counter histories. The difference between the panels is 
most striking. These simulations were designed so that 
all the streams would end with the same shape and or- 
bital phase (although different lengths) in the absence 
of interactions with dark-matter subhalos. Instead, the 



streams show a wide variety of morphologies and orbital 
phases, that can largely be attributed to the few en- 
counters with the more massive dark-matter halos (those 
likely to be observed as stellar satellites). 

Nevertheless, in one way the streams are similar, they 
all show tens of small-scale (~lkpc and below) gaps, 
caused by the many less-massive subhalos — the missing 
satellites. The frequency of these encounters means that 
these disturban ces are apparent in all cases. In this sense, 
we agree with (|Carlbergl 120091 ) that any observed thin 
stream should contain the signatures of missing satel- 
lites. However, it is not clear from Figure [9] that the 
mere existence of old streams rules out the presence of 
a large population of low-mass dark-matter subhalos — 
many examples of thin, old streams survive in our simu- 
lations. Indeed, Figure [9] suggests that we might expect 
to find remnants of destroyed streams at surface bright- 
ness levels below the current level of sensitivity. 

The gaps themselves have a striking morphology: their 
edges are typically not perpendicular to the stream, but 
rather sit at an angle. This "slant" is particularly notice- 
able for debris at the orbital apocenter. We can trace this 
appearance in coordinate space to back to the character 
of the gaps seen in energy /angular-momentum-space in 
earlier figures. The edges of the gaps typically have a sin- 
gle orbital energy with a range of angular-momenta. As 
the debris spreads, the particles sort themselves in energy 
along the stream a nd in angular momen ta perpendicu- 
lar to the stream ([Johnston et al.l[2001h . which results 
in the angular momenta of particles varying monoton- 
ically along the edge of a gap. The slant reflects the 
corresponding monotonic trend in the particles' orbital 
time periods, which are weakly dependent on angular- 
momentum. 

There are also examples of bizarre morphologies in the 
streams in Figure [9] Several streams have discontinu- 
ities as large as several kpc (for example in the middle 
right-hand and bottom left-hand panel) — sufficiently 
disjointed that theses structure might be detected as sep- 
arate streams when we observe them. (Note that the gap 
between the leading and trailing streams in each panel 
is due to the initial conditions in the simulation and not 
subhalo interactions.) 

Some streams even take on a bifurcated appearance, 
for example in the third panel in the middle row and 
the rightmost panel in the bottom row of Figure [9] The 
lower panels of Figure [10] show the final energy / angular- 
momentum distributions for these bifurcated streams, 
with the redder and purpler particles having initially 
higher and lower energy respectively. Comparing this 
to the initial distribution in Figure [2] shows that in both 
cases, the energy-angular momentum distributions have 
been completely flipped by an encounter with a large sub- 
halo, so that the trailing stream ends up on orbits with 
shorter time periods and overtakes the leading stream. 
A bifurcat ion has already been se en in the stellar stream 
from Sgr (Bclokur ov et al.l 120060?) and could possibly be 
due to such an encounter. However, further simulations 
of an Sgr-scale stream (much longer and thicker than 
our Pal-5-like test case) are needed to confirm the plau- 
sibility of this scenario and disting uish it from other ex- 
planations of this bifur cation (e.g. iFellhauer et al.l [20061 : 
iPenarrubia et al1l2010l ) . We do not pursue this investiga- 
tion further here since we are concentrating on encoun- 
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Figure 9. Final particle positions in the X-Y plane for tidal streams evolved in the presence of the full spectrum of ACDM subhalos. 
Each stream begins at differing initial times and positions that were chosen to place all streams at the same ending location in the absence 
other disturbances (see the text). The age is presented in each panel. Note that each panel is centered on the satellite and has a different 
axis scale. 



ters of missing satellites with much colder streams. 

5. DISCUSSION 

Our results suggest that missing satellites, if present 
with number densities predicted by purely dark-matter 
simulations of structure formation in a ACDM Universe, 
should leave their imprint on cold stellar streams such 
as Pal 5 and GD-1 in the form of abrupt changes in sur- 
face density and velocity on few degree-scales and be- 
low. Moreover, this signature should be distinct from 
the few larger-scale deformities expected to be produced 
by known satellites. In this section, we compare our nu- 
merical results to current observations and discuss other 
possible explanations of substructure in stellar streams. 

5.1. The Case of Pal 5: Signatures of Missing Satellites 

The most straight forward thing to observe about stars 
in a stream with current capabilities are their positions 
on the sky. In reality, cold streams like Pal 5 have been 
discovered as overdensities of stars that lie along a re- 
stricted sequence in color-magnitude space (defined by 
a stellar isochrone), rather than direct identification of 
individual stellar members. The identification process 
requires some kind of binning or smoothing of the stars 



in order to decide if a region is overdense, and this will 
soften any i nhomogeneities in the str e am. For example , 
Figure 3 of lOdenkirchen et all (pOOl : iGrillmairl (2009) 
present maps of the region around Pal 5, in which the 
star counts have been smoothed with a parabolic and 
gaussian kernel. For comparison, we select three model 
streams among the ones in Figure [SI count the particles 
in 24' x 24' bins (chosen to mime the level of smoothing in 
the original maps) and present conto urs of surface den- 
sity in Figure HIJ lOdenkirchen et al.1 (|2003l ) found 1,650 
stars in a 10° length of stream, and we expect compara- 
ble surface densities for our particles (of which there are 
3,000 in about 25°). Hence we cho ose the same contour 
levels as lOdenkirchen et al.l (|2003l ) at the 1.5, 2, 3, 5-cr, 
and higher levels, with a — 0.12 particles/arcmin 2 . The 
largest gaps around RA=227°, 243°, and 231° in each 
panel arise from initial conditions (the physical separa- 
tion of the leading and trailing streams) and not subhalo 
encounters. Aside from them, we still find many gaps 
and clumpy substructures in each stream which — while 
somewhat more exaggerated — are in rough agreement 
with the observed Pal 5 stream. 
Moreover, as we discussed in § 14.31 the gaps in the 
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Figure 10. Two streams showing bifurcations arc presented in 
X-Y and energy/angular-momcntum planes. The arrows show the 
direction of the moving streams. The stream particles are color- 
coded by their initial energy (redder higher, purpler lower). The 
crossing points of dashed lines in the bottom panels show the en- 
ergy and angular momentum of initial central particles (marked as 
crosses in the top panels). 



model streams are "slanted" and this morphology is also 
apparent in the contour maps in the middle and right 
panels of Figure II 1 1 (arrows in top and bottom panels in- 
dicated the relevant gaps) . The slants are not as clear in 
the contour maps as in the particle plots since they are 
erased out by smoothing, but can become more promi- 
nent if th e contour l evels are varied. Indeed, in Fig- 
ure 10 of IZou et~a l. (2009), a similar "slanted" gap is 
seen in the bottom-right end of the Pal 5 stream around 
(RA,Dec)=(227°,-2°). The distinctive gap morphology 
caused by subhalo encounters could provide a way of 
distinguishing between competing theories for the ori- 
gin of stream irregularities (see next section). However, 
a higher density of detected stellar traces is needed to 
build a convincing case for the existence of slanted gaps. 

While Figure [TT] indicates that our simulations which 
included the full ACDM mass spectrum appear to be in 
rough agreement with the level of fluctuations observed 
in Pal 5, those that contain only the upper end (i.e. cor- 
responding to known Milky Way satellites) are not so 
successful. Figure [T^] shows the surface density, radial 
velocity dispersion, and radial velocity profiles along the 
streams illustrated in Figure [51 which represent simula- 
tions where the contribution of each mass decade of the 
full mass spectrum was isolated. The profiles arc esti- 
mated from sets each containing 100 particles that have 
been binned in order of increasing angular separation A 
from the central particle. The second and third column 
of panels of Figure [T^] — corresponding to streams dis- 
turbed by low-mass subhalos or "missing satellites" - 
still reveal surface densit y fluctuations which are com- 
parable to the profiles in lOdenkirchen et al.l (|2003[ ) and 



IZou et aLl (|2009h . However, fluctuations of the observed 
amplitudes and scales are not seen in our models that in- 
cluded no satellites (left hand panels), nor in those that 
included only the "visible" satellites from the top few 
decades of the ACDM mass spectrum (right-hand pan- 
els). 

In the future, ideally we would move beyond broad 
consistency with the existence of missing satellites to 
definitively separating and measuring the contribution 
of each decade of the mass spectrum. This would require 
sensitivity to a sufficient number of stars to resolve fluc- 
tuations on sub-degree scales. A comparison of the top 
panels of Figure [T^] with those of Figure [S] illustrates the 
challenge of this with our own simulated data. The dif- 
ference in the nature of fluctuations caused by subhalos 
in the range 10 5 — 10 6 Af Q compared to those caused by 
the mass decade above is striking in the particle plots, 
but the binning of 100 particles to create the profiles 
smoothes over ~ 30° = l°-scalcs so that this signa- 
ture is less apparent in the density profiles. 

Additional information in the form of samples of line- 
of-sight velocity measurements of stars in the region of 
Pal 5's stream are no w being collected, with th e first 17 
already published bv lOdenkirchen et al.l (j2009f ). Figure 
[8j demonstrates that the missing satellites scatter stars 
in velocity space to produce local discontinuities of a few 
km/s on sub-degree scales, while Figure [T2l suggests that 
this scattering may be apparent as a small increase in dis- 
persion in averages over a larger region. However, find- 
ing these signatures observationally remains challenging 
as stream stars cannot be definitively identified and any 
spectroscopic sample will be contaminated by interlop- 
ers. 

5.2. Other Streams 

There are a number of other cold streams with 
small widths and angular extents (and therefore ages) 
that rival or even surpass Pal 5 (for example, GDI, 
NGC5466, Acheron, Cocy t os, Lethe, and Styx streams , 
i Grillmair fc Johnsodl2q06t iGrillmair fc Dionatosll2006bl: 
Fellhauer et al.1 120071 : iGrillmairl 120091 : iKoposov et al.1 
120101) . Detailed maps of these streams' surface densities 
- and of any apparent gaps in them in particular - could 
provide complimentary probes of the missing satellites: 
to first order, since these streams have similar widths and 
all explore the inner few-10 kpc of the Galaxy we expect 
to find similar numbers of and scales to the signatures of 
encounters with dark-matter subhalos. 

Looking ahead, a more ambitious plan would be to 
use multiple streams to examine the relative numbers of 
dark-matter subhalos at different Galactocentric radii. 
This interpretation would require rather better than sim- 
ple intuitive extrapolations from the current study for the 
number a nd scales of gaps. For example, from the esti- 
mates in IKoposov et all (|2010f ). while both Pal 5 and 
GDI debris sits at similar Galactocentric radii, GDI 
orbits slightly farther out (between 18-30kpc from the 
Galactic Center) and should perhaps therefore encounter 
fewer missing satellites and contain fewer gaps. However, 
although GDI has a similar width to Pal 5, it has a much 
longer angular extent (~ 70°, or ~ 40° as viewed from 
the Galactic center): naively, a longer extent suggests 
GD 1 is the older stream, which points in the direction 
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Figure 11. Particle distribution (upper panels) and surface density maps (lower panels) for three model streams from Figure [9] projected 
onto the plane of the sky. In the particle plots, many slanted gaps are found which are in some cases still apparent in the contour maps 
(indicated by arrows). 
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Figure 12. Top-Bottom: Number density, velocity dispersion, and radial velocity profiles along the stream. Left-Right: The stream 
evolved without subhalos, and with subhalos in the mass range 10 5 - 10 6 /10 6 - 10 7 /l(y - 10 s /10 8 - 1O 9 M are presented. A is the angular 
separation from the host satellite center. Negative and positive A mean trailing and leading stream respectively. Mean errors in the counts 
in the upper panels are indicated by the bar in the top-right side of the panels in the top row. 
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of more gaps. To further complicate things, GDI is at 
pericenter, while Pal 5 is at apocenter, and gap scales 
and morphologies are phase-dependent. All these differ- 
ences would have to be taken into account with specific 
models of each stream before any meaningful comparison 
could be made. 

5.3. Other Influences on Stream Non-uniformity 

There are several other processes that would affect the 
appearance of clumps in stellar streams which were not 
included in our numerical tests. In some cases, these 
could plausibly provide alternative explanations of clump 
origins which do not require the existence of missing 
satellites. 

: Multiple or continuous mass-loss events — In our sim- 
ulations, the test particles all started at the same 
time and orbital phase, mimicking the evolution of 
debris form a single mass-loss event and resulting 
in debris in which the orbital properties of particles 
varied monotonically with angle along the stream. 
As a consequence, subhalos in our simulations di- 
rectly impacted particles that were not only close 
physically, but also close in their orbital properties. 
In reality, streams are likely to be composed of stars 
lost over extended periods and during several mass- 
loss events, so that subhalos would affect debris 
with a wider range of orbital properties. Hence the 
gaps in energy space created by encounters would 
not be as clean as those seen in our figures (for 
example in Figure [5]) . This in turn suggests that 
our simulations have likely overestimated the am- 
plitude (though not physical scales) of fluctuations 
in the final spatial distribution of the debris: more 
realistically, stars from different mass-loss popula- 
tions would overtake each other across the spatial 
gaps and soften the contrast. 

: Asymmetries in the global potential — We assumed 
a purely spherical parent potential, neglecting as- 
phcrical contributions from the disk, bar, spiral 
arms, and Galactic warp or the (likely) triaxiality 
of the dark-matter halo. iSiegal-Gaskins fc Valluril 
(2008) have shown in particular how the last of 
these can distort (relatively hot) stellar streams 
from large Galactic satellites. However, we antic- 
ipate that none of these large asymmetries would 
leave signatures that mimic the degree- and sub- 
degree scale gaps seen in our experiments on colder 
stellar streams and specifically due to direct hits by 
missing satellites. 

: Other sources of gravitational shocks — Streams can 
also suffer shocks from pericentric encounters with 
the Galactic bar/bulge and when passing through 
the Galactic disk. Disk shock, which are not taken 
into acc ount in our models, can destroys globular 
clusters (|Gnedin et al.lll999l) . but are not believed 
to be the source of discontinuities in Pal 5's stella r 
streams ()Dehnen et al.l 12004 IKupper et al.ll2010[ ). 
The effect of a disk passage could be enhanced if 
the stream happened to pass directly through a Gi- 
ant Molecular Cloud (GMC). By extrapolating the 
GMC mass-function, we estimate there to be of or- 
der one thousand of CMC's within the Solar Circle 



in the disk with masses < 10 7 M ( t 1 (|Solomon et al.l 
119871: iWilliams fc McKe"elll997D . which could imply 
a greater density in this region than dark-matter 
subhalos of comparable scales. However , although 
Pal 5 crosses the Galactic plane 35 times during the 
8.44 Gyrs of our simulation, it only passes through 
the inner disk(< 8kpc) 8 times and the closest ap- 
proach to the Galactic center when it crosses the 
disk is 7.5kpc. Hence, far less encounters of GMCs 
than the number of direct hits from missing satel- 
lites are expected. Overall, we conclude that the 
gravitational shocks from these sources should be 
much less frequent than encounters with subhalos 
for the Pal 5 stream. 

: Self gravity of the stream — Our test-particle exper- 
iments do not capture the Jean's instability that 
could cause clumpy substr ucture in a stream even 
in th e absence of subhalosf Q uillen fc Comparettal 
I2010D . The importance of this effect depends on the 
density and velocity dispersion of the stream as well 
as the eccentricity of the stream orbit. Simulations 
including stream self-gravity are needed to assess 
whether it is possible to disentangle the signatures 
of subhalo encounters from Jeans instability. 

: Epic yclic m o tions of es c aping s tars — IKupper et all 
(|2008l ): Uust et al.l (f2009T) ; IKupper et all (|2010f ) 
used analytic estimates and numerical experiments 
to demonstrate that fluctuations in stream density 
can arise from epicyclic motions of escaping stars. 
Their experiments showed that these fluctuations 
should be most apparent in cases where the stars 
escape "evaporatively" (at low-velocity relative to 
the satellite) and the debris is observed at pericen- 
ter. At apocenter of their orbit with eccentricity 0.5 
(comparable to Pal 5), the fluctuations appear to 
be less than a factor of two in amplitude and rather 
regularly spaced, in clear contrast to the prominent 
irregularly-spaced gaps induced in streams by miss- 
ing satellites in our experiments. 

Overall, we anticipate that the mixing of debris from 
multiple mass loss events could obscure the contrast and 
distinct morphology of the gaps seen in our idealized sim- 
ulations. However, these softened fluctuations might ac- 
tually be more consistent with those seen in observations 
of Pal 5's stream. Of the remaining missing effects, the 
self-gravity of the stream seems to be the dominant factor 
confusing our interpretation and weakening our conclu- 
sion that the "missing satellites" have been found. 

6. CONCLUSIONS 

The goal of this study was to test if a dynamically 
cold stellar stream could survive in the presence of dark- 
matter subhalos and if so, to characterize observable fea- 
tures in the stream produced by subhalo encounters. In 
particular we were interested in finding signatures of sub- 
halos that would be considered "missing satellites" rather 
than those that we already know exist. We conclude that: 

1. The mere existence of cold stellar streams does not 
imply the absence of missing satellites — dynam- 
ically cold streams can survive for many Gyears 
even when bombarded by subhalos. 
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2. Those streams that are observed should contain the 
imprint of past direct impacts from subhalos in the 
form of gaps in surface density and discontinuities 
in velocities. 

3. The frequency and scale of the gaps is dependent 
on the mass spectrum of subhalos and the prop- 
erties of the stream itself. In the case of Pal 
5, there should be observable fluctuations at de- 
gree and sub-degree scales due to ~ 90 encoun- 
ters with subhalos in the mass range 10 5 — 10 7 M Q 
(i.e. the missing satellites), while distant encoun- 
ters with larger subhalos produce less observable 
effects. Hotter stellar streams, such as Sgr's debris 
are large enough to hide the signatures of the many 
encounters it suffers with missing satellites. 

4. Current observations of Pal 5's stream show that 
its surface density profile contains fluctuations on 
comparable scales to those predicted in our simu- 
lations, and could be interpreted as direct proof of 
the existence of missing satellites. 

Our study has built on previous works to look at: ef- 
fects of individual encounters with different mass scales; 
the distinct signatures of subhalos in different mass 
decades; the integrated influences of subhalo encounters; 
and comparisons to the current observational data. How- 
ever, while we have succeeded in outlining the expected 
scales of signatures of different decades of the ACDM 
mass spectrum in stellar streams, there are still some 
open questions that need to be answered before we are 
confident enough to say that we have definitely found 
the missing satellites. In particular, self-consistent sim- 
ulations of Pal 5 disruption including stream self-gravity 
and continuous mass loss are needed to clarify to what 
extent subhalo signatures could be confused by these ef- 
fects. With such simulations in hand, it would then make 
sense to develop a multi-dimensional statistic to compare 
to all available data in a more robust way. The promise 
of current and near-future efforts to map Pal 5 and other 
stellar streams are strong motivators for this work. 
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